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ABSTRACT 

Precise spectra of 3C 279 in the 0.5-70 keV range, obtained during two epochs of 
Swift and NuSTAR observations, are analyzed using a near-equipartition model. We 
apply a one-zone leptonic model with a three-parameter log-parabola electron energy 
distribution (EED) to fit the Swift and NuSTAR X-ray data, as well as simultane¬ 
ous optical and Eermi-LAT y-ray data. The Markov Chain Monte Carlo (MCMC) 
technique is used to search the high-dimensional parameter space and evaluate the 
uncertainties on model parameters. We show that the two spectra can be successfully 
fit in near-equipartition conditions, defined by the ratio of the energy density of rela¬ 
tivistic electrons to magnetic field Ce being close to unity. In both spectra, the observed 
X-rays are dominated by synchrotron-self Compton photons, and the observed 7 rays 
are dominated by Compton scattering of external infrared photons from a surrounding 
dusty torus. Model parameters are well constrained. From the low state to the high 
state, both the curvature of the log-parabola width parameter and the synchrotron 
peak frequency significantly increase. The derived magnetic fields in the two states 
are nearly identical 1 G), but the Doppler factor in the high state is larger than 
that in the low state (~28 versus ^18). We derive that the gamma-ray emission site 
takes place outside the broad-line region, at > 0.1 pc from the black hole, but within 
the dusty torus. Implications for 3C 279 as a source of high-energy cosmic-rays are 
discussed. 

Key words: galaxies: jets — gamma rays: galaxies — radiation mechanisms: non- 
thermal 


1 INTRODUCTION 

Blazar emission is generally dominated by non-thermal ra¬ 
diation over all frequencies ranging from radio to TeV 7 
rays. By modeling blazar spectra, emission mechanisms and 
physical properties of the relativistic jets can be determined. 
The typical multi-wavelength spectral energy distribution 
(SED) of a blazar is characterized by two distinct humps. 
It is generally accepted that the first hump in the blazar 
SED is nonthermal synchrotron emission radiated by rela¬ 
tivistic electrons in the jet. The origin of the emission in 
the 7 -ray hump is less certain. Leptonic and hadronic mod¬ 
els have both been proposed to explain the formation of 
the second bump (see, e.g., Bottcher et al. 2012, for review). 
In general, both types of models are able to reproduce the 
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SEDs well, but require quite different jet properties, with the 
hadronic models sometimes requiring super-Eddington jet 
powers (Bottcher et al. 2013; Zdziarski & Bottcher 2015). 
The leptonic models are more attractive by allowing much 
smaller jet powers, and by more readily explaining rapid 
blazar variability from highly radiative electrons in con¬ 
trast to the more weakly radiating hadrons (cf. Cerruti et al. 
2015). 

All hadronic models are hybrid lepto-hadronic, because 
the synchrotron bump is believed to be radiated by di¬ 
rectly accelerated electron/lepton primaries in the jet. Be¬ 
sides the leptonic synchrotron self-Compton (SSC) compo¬ 
nent, 7 rays in hadronic models can result from proton- 
or ion-synchrotron radiation (Aharonian 2000; Miicke et al. 
2003) and p 7 interactions, including both photopion 
and Bethe-Heitler photopair production (Atoyan & Dermer 
2001, 2003; Bottcher et al. 2009, 2013; Murase et al. 
2014; Mastichiadis et al. 2013; Dimitrakoudis et al. 2014; 
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Weidinger & Spanier 2015; Yan & Zhang 2015). Interac¬ 
tions of ultra-relativistic protons and ions in the jet with 
photons indnce cascades by the injection of extremely high- 
energy photons and leptons. 

In leptonic models, 7 rays are prodnced via in¬ 
verse Compton (IC) scattering of the relativistic elec¬ 
trons, including SSC scattering (e.g., Maraschi et al. 
1992; Tavecchio et al. 1998; Finke et al. 2008; Yan et al. 
2014) and external Compton (EC) scattering (e.g., 
Dermer & Schlickeiser 1993; Sikora et al. 1994; 
Blazejowski et al. 2000; Dermer et al. 2009; Kang et al. 
2014; Paliya et al. 2015). EC models seem to be required 
for flat-spectrum radio quasars (FSRQs; a subclass of 
blazars with intense broad-line radiation in the inner-jet). 
Depending on the emission site, the external seed photons 
could be dominated by radiation from accretion-disk 
(e.g., Dermer et al. 1992; Dermer & Schlickeiser 1993, 
2002), broad-line region (BLR) (e.g., Sikora et al. 1994; 
Zhang et al. 2012), and dusty torus (e.g., Blatejowski et al. 
2000). Many efforts have been made to locate the 7 -ray 
emission regions of FSRQs (e.g., Liu & Bai 2006; Bai et al. 
2009; Tavecchio et al. 2010; Yan et al. 2012; Dotson et al. 
2012; Nalewajko et al. 2014), but it is still an unresolved 
issue. 

Radio telescopes can partially resolve the structure of 
blazar jets. The observed connections between 7 -ray flares 
and radio/millimeter flares are used to argue for gamma-ray 
emission arising at distance scales of > pc, or even > 10 pc 
(e.g., Jorstad et al. 2001, 2010; Agudo et al. 2011). However, 
Tavecchio et al. (2010) argued that the rapid variabilities on 
timescale of a few hours challenge the scenario in which 7 -ray 
emission site is placed at tens of parsec from the black hole. 
Using the time lags of 7 rays relative to broad emission lines, 
Liu, Bai & wang (2011) argued that 7 -ray emission may be 
produced in the regions of the jet at < 2 pc from the black 
hole. 

Although blazar emission models have been constructed 
that accurately reproduce contemporaneous SED data, the 
number of free parameters in even the simplest SSC model 
can exceed the number of constraining observables. There¬ 
fore, blazar SED modeling generally yields degenerate pa¬ 
rameter values, which makes it difficult to investigate corre¬ 
lations between model parameters. 

Claims and explanations of correlations between model 
parameters have been made, for example, for statistical 
quantities like the blazar sequence (Fossati et al. 1998; 
Ghisellini et al. 1998), or for flux-flux or flux-index corre¬ 
lations between different spectral states of the same blazar. 
Due to many free and unconstrained parameters, it is diffi¬ 
cult to establish robust correlations between model param¬ 
eters using simple fitting methods. 

Using the Markov Chain Monte Carlo (MCMC) tech¬ 
nique, Yan et al. (2013) constrained parameter space in an 
SSC model for Mrk 421. The high-quality SEDs enable us 
to derive the confidence intervals of each model parameter 
using the MCMC method, and to do a comparison of model 
fits for different SEDs. In a series of papers (Yan et al. 2013; 
Peng et al. 2014; Zhou et al. 2014), we showed that fitting 
high-quality SEDs with the MCMC technique is a powerful 
approach to investigating the blazar jet physics. In a related 
modeling effort, Dermer et al. (2014) assumed that a three- 
parameter log-parabola function approximates the electron 


energy distribution (EED). Introducing the equipartition 
parameter = u^^jug, as the ratio of the energy den¬ 
sities of electrons and magnetic field VLg, gives, when 
((e « 1, the minimum synchrotron jet power. Besides ^e, the 
other three equipartition factors in this modeling approach 
(Cerruti et al. 2013; Dermer et al. 2014) are: (1) Cs, the ra¬ 
tio of the energy density of synchrotron radiation to 
v!gi\ (2) i(blr, the ratio of the energy density of BLR radia¬ 
tion in comoving frame Wblr lo and (3) ((dust, the ratio 
of the energy density of dust radiation in comoving frame 
Wdust to u'gi. In this approach, the physical parameters (e.g., 
fluid magnetic field B' and Doppler factor 5 d) are expressed 
in terms of observables (e.g., synchrotron peak frequency 
and synchrotron peak luminosity, and variability time scale). 
With this formulation, GeV cutoffs in blazar spectra are 
explained under equipartition conditions by IG scattering 
broad-line region (BLR) radiation in Klein-Nihisna regime 
(Cerruti et al. 2013). In SSC models, the equipartition mod¬ 
eling approach allows the complete set of parameters to be 
fully constrained, contingent on the underlying assumption 
that the electron distribution can be described, at least ap¬ 
proximately, by a log-parabola function. Here, we constrain 
the parameters using the near-equipartition, log-parabola 
model for 3C 279. In our approach, we do not adopt ((blr 
and ((dust to define the energy densities of the external radi¬ 
ations. Instead, we use the distance r from emitting region 
to super-massive black hole to deduce the energy densities 
of the external radiations via the standard scaling relations 
for BLR and dust torus radiations. We will give the details 
in Section 2. 

The blazar 3C 279, at redshift 2 = 0.536, is classified 
as an FSRQ due to its strong prominent emission lines, and 
was one of the first 7 -ray blazars discovered (Hartman et al. 
1992). Its synchrotron peak frequency is < Hz, mak¬ 
ing it a low-synchrotron peaked object, like most FSRQs 
(Ackermann et al. 2011). Rapid variability is observed in 3C 
279, the analysis in Paliya (2015) shows that the variability 
timescale for GeV emission of 3C 279 can be down to - 
2 hour. 

Using the near-equipartition leptonic model with a log- 
parabola electron energy distribution (EED), Dermer et al. 
(2014) satisfactorily explained four SEDs of 3C 279 reported 
in Hayashida et al. (2012). Since then, detailed simultaneous 
multiwavelength SED data of 3C 279 in 4 states of elevated 
activity that took place between 2013 December to 2014 
April were analyzed by Hayashida et al. (2015). In particu¬ 
lar, Hayashida et al. (2015) reported 0.5-70 keV X-ray spec¬ 
tra obtained from Swift-XKT and NuSTAR observations in 
two flaring periods, namely Period A, taking place from 16 
- 19 December 2013, and Period C, from 31 December - 3 
January 2014. Hayashida et al. (2015) argue that their one- 
zone leptonic models with a double-broken power-law EED 
fail to explain the Stui/t-XRT and NuSTAR X-ray spectra. 

In this paper, we fit the two simultaneous SEDs that in¬ 
clude NuSTAR data reported in Hayashida et al. (2015) us¬ 
ing the near-equipartition blazar model with a log-parabola 
EED. Using the MCMG technique, we derive the best-fit re¬ 
sults and the uncertainties on parameters. The implications 
of the results on the acceleration processes in the jet of 3C 
279 are discussed. Supposing that 3C 279 is a high-energy 
cosmic-ray (CR) source, we discuss its implied CR luminos- 


MNRAS 000, 1-10 (2015) 


Parameter constraints in Modeling of SC 279 3 


ity. We use parameters Hq — 71 km s ^ Mpc = 0.27, 

and r^A = 0.73. 


2 EMISSION MODEL AND FITTING 
TECHNIQUE 


In the log-parabola leptonic blazar model (Dermer et al. 
2014), the non-thermal electron distribution is assumed to 
be isotropic in the emission region (blob), and described by 
a log-parabola function 





log(7'/7pk) 


( 1 ) 


where p is the electron Lorentz factor, b is the spectral cur¬ 
vature parameter, and 7 pk is the peak Lorentz factor in the 
distribution. Emission from the blob is strongly 
boosted due to the beaming effect. Besides SSC emission, EC 
components for low-energy target photons from the broad¬ 
line region (BLR; EC-BLR) and the dusty torus (EC-dust) 
are included. 

As mentioned in Section 1, we do not use the 
equipartition relations ((blr and ^dust (Cerruti et al. 2013; 
Dermer et al. 2014) to normalize the energy densities of the 
external radiations. We take the distance r from the central 
black hole to the emitting blob as an input parameter. The 
energy densities of BLR radiation (ublr) and dust radiation 
(udust) can be expressed as functions of r (Sikora et al. 2009; 
Hayashida et al. 2012) 


ItBLR(r) 

TBLRTdisk 

(2) 

47rr|,LRc[l -h (r/rBLR)3] ’ 

Udust (yC ) 


(3) 

47rr-i„gtC[l-b (r/rdust)^l ■ 


The size of BLR is related to the disk luminos¬ 
ity Ldisk: rBLR = 10^’'(Ldisk/10'‘® erg s“^)^/^ cm 

(Ghisellini & Tavecchio 2009; Ghisellini et al. 2014). We 
assume a dust torus with the size of rdust = 
10^®(Ldisk/10^® erg s“^)^^^ cm. This size of dust torus is half 
of the size used in Ghisellini et al. (2014), which enhances 
the energy density of dust radiation. Using the scalings be¬ 
tween I/disk and rBLR/dust, we can rewrite Eqs. (2) and (3) 
as 


ublrO 


'^dust (c) 


0.3tblr 
1 -I- (r/rBLR)^ 


0.003rdust 
1 + (r/rdust)'‘ 


-3 

erg cm , 

(4) 

erg cm“®. 

(5) 


where tblr and rdust are the fractions of the disk luminosity 
reprocessed into BLR radiation and into dust radiation, re¬ 
spectively. We adopt the typical values of tblr = 0.1 (e.g., 
Ghisellini et al. 2014) and rdust = 0.3 (e.g., Hao et al. 2005; 
Malmrose et al. 2011). Using Eqs. (4) and (5), one can ob¬ 
tain UBLR(erg cm“®) < 0.3rBLR, and Udust(erg cm“®) < 
Srdust X 10“®. The two energy densities are limited by tbur 
and rdust, respectively. When r is smaller than rdust, Udust 
varies by a factor of at most 2. Moreover, the energy density 
of the IR dust radiation should be lower than the energy den¬ 
sity of a blackbody with the temperature of Tdust, namely 
Mdust < Wbb(2dust) 3 X 10“'*(rdust/440 K)'‘. 


The external radiation is assumed to be a dilute black- 
body radiation, namely a blackbody spectral shape normal¬ 
ized to UBLR/dust- BLR radiation is dominated by Lya line 
photons. We adopt an effective temperature for the BLR 
radiation of Tblr = 6.3 x 10“^ K, so that the energy den¬ 
sity of BLR radiation peaks at « 2.S2kBT/h = 3.7 x 10^® 
Hz in the u{v) distribution (Tavecchio & Ghisellini 2008). 
We consider IR dust radiation with Tdust = 800 K. Note 
that Dermer et al. (2014) approximated the BLR and dust 
photon fields as monochromatic. Cerruti et al. (2013) con¬ 
sidered a complex of atomic emission lines in BLR. The dif¬ 
ferent approximations for the BLR and dust radiations do 
not significantly modify the modeling results. 

The input parameters in the model are; (i) — 

Lsyn/lO"*® erg s“^, the apparent isotropic bolometric syn¬ 
chrotron luminosity; (ii) zzi 4 = {1-\- /iC'^ Hz, the syn¬ 

chrotron peak frequency in the source frame, where 
is the measured synchrotron peak frequency; (iii) = 
tobs.var/[(l-l- 2 ) 10 '* s)], the source variability timescale, where 
tobs.var is the minimum measured variability timescale; as 
already noted, we define the equipartition parameter (iv) 
Ce = Ug/us; (v) Cs = Usy^/u's, the ratio between the syn¬ 
chrotron photon (w^yn) and magnetic-field energy densities; 
(vi) b, curvature parameters of EED; (vii) r, the location of 
emitting blob along the jet; (viii) Tdust, the temperature of 
dust torus IR radiation; and (ix) Ldisk, accretion disk lumi¬ 
nosity. 

There are six output parameters: (1) Doppler factor 
5 d; (2) fluid magnetic-field strength B'- (3) (4) ublr; 

(5) Udust; and (6) comoving radius of emitting blob, R' = 
c5Dtobs,var/(l + z). The physical model parameters, 5 d, B', 
and 7pk, are deduced by using the equations in Dermer et al. 
(2014, 2015). The SSC and EC emissions are calculated us¬ 
ing the methods given in Dermer et al. (2009). Synchrotron 
self-absorption (SSA) is included. The contribution of ther¬ 
mal emission is included, which is assumed to be from a 
Shakura-Sunyaev disk, using the expression following eq. (9) 
in Dermer et al. (2014). 

We use the MCMC technique to do the fitting (see de¬ 
tails in Yuan et al. 2011; Liu et al. 2012; Yan et al. 2013). 
The MCMC method, based on the Bayesian statistics, is a 
powerful tool for high-dimensional parameter space investi¬ 
gation. The Metropolis-Hastings sampling algorithm which 
ensures that the probability density functions of model pa¬ 
rameters can be asymptotically approached with the num¬ 
ber density of samples, is adopted to determine the jump 
probability from one point to the next in parameter space 
(Gamerman 1997; Mackay 2003). The MCMC method is also 
an effective approach for determining uncertainties of model 
parameters. 


3 RESULTS 

We consider the SEDs of 3C 279 from simultaneous obser¬ 
vations of Swift, NuSTAR, and Termi-LAT in two periods 
reported in Hayashida et al. (2015), namely Period A (2013 
December 16-19) and Period C (2013 December 31-2014 Jan¬ 
uary 3). Since an X-ray spectrum is lacking in Period B and 
there is no NuSTAR data in Period D, we do not consider 
these SEDs. Note that the X-ray flux increased by a factor 
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Figure 1. Best-fit model to the SED of 3C 279 for Period A 
(data from Hayashida et al. 2015). The inset shows the details of 
fit at X-ray energies. 

of « 2 and the 7 -ray flux increased by a factor of « 3 from 
Period A to Period C. 

The X-ray flux showed no intraday variations in 
Period A, but showed intraday variations in period C 
(Hayashida et al. 2015). Therefore, the variability timescale 
for Period A is taken to be t 4 = 5 which corresponds to an 
observed variability timescale of ^^21 hours, and the vari¬ 
ability timescale for Period C is set to be f 4 = 1 which cor¬ 
responds to an observed variability timescale of ~4 hours. 

There is no evidence for a thermal emission feature in 
the two SEDs, making it difficult to constrain the accretion 
disk luminosity Ldisk- We adopt Ldisk = 1.5 x ergs“^, 
which is the maximum disk luminosity allowed by the 
optical-UV SED. Using this value, we derive ruLR — 0.04 pc 
and rdust ^ 0.4 pc. 

By comparison, at a much earlier 3-week epoch be¬ 
tween 1992 December and 1993 January, Pian et al. (1999) 
originally discovered the accretion disk-emission from 3C 
279 during a low-state monitored with ROSAT, lUE, and 
EGRET, finding a uF^ flux of « 3 x 10“^^ erg s“^ cm“^ be¬ 
tween « 1 and 2x 10^® Hz, corresponding to a UV luminosity 
Ddisk = 2.4 X 10^® ergs"^ 

Due to the strong SSA below 10^^ Hz, we fit only the 
optical. X-ray and 7 -ray data, but not the 230 GHz SMA 
point. 

Fig. 1 shows the best-fit to the SED of 3C 279 for Period 
A. The fit is satisfactory. The joint Swift and NuSTAR X- 
ray spectrum is successfully fit, and is dominated by the 
SSG component. The 7 -ray emission is dominated by the 
EG-dust component. 

In Fig. 2, we give the one-dimensional probability dis¬ 
tributions and the two-dimensional contours of input and 
output parameters. The complete information on constrain¬ 
ing parameters can be read from this figure. We note that 
(^e is constrained in the range of [0.8 - 1.7] (the marginalized 
95% confidence interval). The distance r is constrained in the 
range of [0.1 - 0.4] pc, namely the emitting blob is outside 


the BLR, but is still inside the dust torus (rsLR — 0.04 pc 
and rdust — 0.4 pc). The derived Udust is in the range of 
[5 - 9[xl0“'* erg cm“®, and ublr is in the range of [0.2 - 
10]xl0“‘^ erg cm“®. The magnetic field B' is in the range of 
[0.7 - 1.2] G, and Jd is in the range of [17 - 20]. We summarise 
the marginalized 95% confidence intervals for input and out¬ 
put parameters in Table 1. We also notice that there is a cor¬ 
relation between r and C,e\ and B' is anti-correlated with Jd 
( see the two-dimensional contours). These parameter cor¬ 
relations showed in Fig. 2 are caused by degeneracies of 
the leptonic model (e.g., Tavecchio et al. 1998; Sikora et al. 
2009; Dermer et al. 2014). Nevertheless, all parameters ex¬ 
cept ublr are well constrained. 

The fit to the SED from Period C is also satisfactory 
(Fig. 3). The Swift-NuSTAR X-ray spectrum from Period C 
shows a break at 3.7 keV with photon indices of 1.37 and 
1.76, respectively, below and above the break energy, when 
fit by a broken power law (Hayashida et al. 2015). Our re¬ 
sults show that the one-zone leptonic model can successfully 
fit the X-ray spectrum. The X-ray emission in Period G is 
almost entirely attributed to SSG, and the spectral break is 
naturally explained. The 7 -ray emission is again dominated 
by the EG-dust component. 

The one-dimensional probability distributions and the 
two-dimensional contours of input and output parameters 
derived in the fitting to the SED for Period G are shown 
in Fig. 4. In this Period, the 7 -ray emission site r is in the 
range of [0.1 - 0.4[ pc, identical to that in Period A; but a 
larger of [3 - 5] and Jd of [26 - 30[ is required to explain 
the higher 7 -ray flux. The magnetic field B' of [0.9 - 1.4] 
is nearly identical to that in Period A. From Period A to 
Period G, and vn increase, respectively, from -^0.5 to 0.8, 
and from -^0.3 to 0.5; and b increases from ~1.2 to 1.5. 

In Table 1, we also give the jet powers for a two-sided jet 
in the form of radiation (Pr), Poynting flux (Pb), where we 
assume that Doppler factor is related to the bulk Lorentz 
factor through Jd = P. The relativistic emitting electrons 
power Pe, which is not shown, is obtained from the relation 
Pe=^ePB. This is the minimum jet power, and the addition 
of hadrons will only increase this power. One can see that 
the radiation power Pr is significantly greater than Pb and 
Pe in both cases. 


4 DISCUSSION AND CONCLUSIONS 

A one-zone leptonic model with a three-parameter log- 
parabola EED (Dermer et al. 2014) is used to fit the two 
SEDs of 3C 279 where joint NuSTAR and Swift data 
are available (the SEDs from Period A and Period G in 
Hayashida et al. 2015). Using the MGMC method, we ob¬ 
tained the best-fit results and the uncertainties on the input 
and output parameters. Our results show that the two SEDs 
can be successfully fit in near-equipartition conditions. In 
both cases, the 7 -ray emission is dominated by the EG-dust 
component, and the EC-BLR component is essentially negli¬ 
gible; and the X-ray emission is dominated by SSG emission. 
The SSG origin of X-ray emission is expected to be highly 
polarized if the magnetic field is perpendicular to the line of 
sight (e.g., Krawczynski 2012). The X-ray measurements of 
the degree of polarization by future detectors such as X-ray 
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Figure 2. One-dimensional probability distribution of parameter values (left; the dashed curves are the mean likelihoods of samples and 
the solid curves are the marginalized probabilities) and two-dimensional contours of parameters (right; the regions enclosing 68%(95%) 
confidence level are shown) for Period A. The upper panel is for the input parameters, and the lower panel is for the output parameters. 


Table 1. Input and output parameters values. The mean values and the marginalized 95% confidence intervals (Cl) for interested 
parameters are reported. 


Input 


Ce 

b 

-^48 

1/14 

*4 

Cs 

r 

^dust 









(pc) 

(K) 

(10'^*^ erg s -*•) 

Period A 

1.12 

1.24 

0.18 

0.30 

5 

0.47 

0.2 

800 

0.15 

(95% Cl) 

0.80-1.71 

1.15-1.32 

0.17-0.19 

0.26-0.34 

- 

0.41-0.52 

0.1-0.4 


- 

Period C 

4.06 

1.49 

0.23 

0.54 

1 

0.81 

0.3 

800 

0.15 

(95% Cl) 

3.04-5.0 

1.43-1.57 

0.22-0.24 

0.50-0.58 

' 

0.74-0.89 

0.1-0.4 


- 

Output 


B 

■^D 


'^dust 

■^BLR 

r' 

P-B 

Pt 



(G) 



(10~^ erg cm~^) 

(10~^ erg cm~^) 

(10^® cm) 

(10^® erg s-1) 

(lO"^® erg s“^) 


Period A 

1.0 

18 

260 

0.7 

0.2 

2.7 

1.7 

2.0 


(95% Cl) 

0.7-1.2 

17-20 

220-290 

0.5- 0.9 

0.02-1 

- 

- 

- 


Period C 

1.1 

28 

300 

0.6 

0.1 

0.8 

0.5 

1.3 


(95% Cl) 

0.9-1.4 

26-30 

280-310 

0.5- 0.9 

0.03-2 

- 

- 

- 



timing and polarization (XTP) and ASTRO-H would give 
more details on the magnetic field. 

We take advantage of the MCMC technique to search 
the multi-dimensional parameter space, and to evaluate the 
uncertainties on the parameters. We have shown that all pa¬ 
rameters except ublr are constrained very well by the cur¬ 


rent data. The stringently constrained parameters enable us 
to investigate the important issues in blazar physics confi¬ 
dently. 

We derive the energy density of BLR radiation ublr < 
0.002 erg cm~® at 95% confidence level. Our results show 
that the EC-BLR component is essentially negligible in mod- 
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Figure 3. Best-fit model to the SED of 3C 279 for Period C (data 
from Hayashida et al. 2015). The inset shows the details of fit at 
X-ray energies. 


elling the two SEDs 3C 279. Note that Cerruti et al. (2013) 
showed that EC-BLR component is necessary to explain the 
GeV break in 3C 454.3. The distance r from 7 -ray emission 
blob to the black hole is well constrained in the range of [ 0.1 - 
0.4] pc, namely the 7 -ray emission region is outside the BLR, 
but is inside the dusty torus. The value of r depends on the 
assumptions on recR and rdust. Compared with the previous 
works (e.g., Ghisellini & Tavecchio 2009; Hayashida et al. 
2012; Nalewajko et al. 2014; Ghisellini et al. 2014), we 
adopted a smaller ruust- In Appendix A, we show the fit¬ 
ting results with ruust = 2 x 10 ^®(I/disk/ 10 ^® erg s“^)^^^ cm 
(Ghisellini et al. 2014). We find that the lower limit of r 
derived in the fittings with the larger rdust is still 0.1 pc; the 
upper limit is modified, but is still less than rdust. Therefore, 
our conclusion on the 7 -ray emission site is still tenable. Our 
result is consistent with the 7 -ray emission site in 3C 279 
derived by Dermer et al. (2014), and is also consistent with 
the result derived by Nalewajko et al. (2014) who used an 
independent method to constrain the 7 -ray emission site. 
Pacciani et al. (2014) explored the emission zone in high- 
energy flares of 10 FSRQs, and also found the evidence of 
gamma-ray flares occurring outside the BLR. 

The log-parabola EED can be generated by stochas¬ 
tic particle acceleration, for example, through systematic 
gyro-resonant acceleration of electrons with plasma waves 
(e.g., Becker et al. 2006; Tramacere et al. 2011; Asano et al. 
2014; Kakuwa et al. 2015), and the correlations between the 
curvature b and the peak energy 7 pij of EED provide evi¬ 
dence about the acceleration process (e.g., Tramacere et al. 
2011). By fitting the SEDs, the EEDs are inferred, which 
helps to identify acceleration processes in the jet (Yan et al. 
2013; Peng et al. 2014; Zhou et al. 2014). Indeed, Yan et al. 
(2013) found that, because of the improved fits with curv¬ 
ing rather than power-law EEDs, it is likely that stochas¬ 
tic acceleration is acting in the jet of Mrk 421 in the gi¬ 
ant flare state in February 2010. The acceleration process 
is also revealed by dynamic changes of the SED, which im¬ 


plies changes in the EED. Yan et al. (2013) also found that 
for Mrk 421 as V 14 increased from ~ 1000 in the low state 
to ~ 10000 in the giant flare, the curvature of PLLP EED 
in N'{y) distribution increased from ~ 1.7 to ~ 3.8. We 
discussed that the changes in and the curvature in Mrk 
421 are not compatible with a purely acceleration-dominated 
scenario. The two states of 3C 279 that we have analyzed 
indicates that a changing EED in terms of a changing log- 
parameter width parameter b plays an important role in dif¬ 
ferent spectral states. 

It is known that X-ray spectra of several high- 
synchrotron-peaked (HSP) BL Lac objects can be described 
by a log-parabola function (e.g., Massaro et al. 2004a,b). 
Using a large data set of X-ray observations for several 
HSPs with > 10, Tramacere et al. (2007, 2009) and 
Massaro et al. (2008) find that the synchrotron peak fre¬ 
quency is anti-correlated with the curvature parameter ob¬ 
tained by fitting the X-ray spectrum. Chen (2014) finds that 
the synchrotron peak frequency is anti-correlated with the 
curvature parameter of the synchrotron spectrum for a sam¬ 
ple of Fermi blazar. 

Tramacere et al. (2011) discussed the vu - b trends in 
an SSC model with stochastic acceleration. They showed 
that the anti-correlation between i>i 4 and b for HSPs could 
be explained in the stochastic acceleration SSC model by a 
change in the diffusion process rather than by a change of 
magnetic field that affects the cooling. The radiative cooling 
of electrons in FSRQs is more efficient than that in HSPs, 
which may lead to the different 1/14 - b trend from that 
for HSPs. Our results show that both 1/14 and b significantly 
increase from period A to period C. More high-quality SEDs 
are needed to confirm the V 14 ~ b trend for FSRQs. Moreover, 
our results indicate that the 7 -ray emission region is in the 
dust torus, in contrast to Zhang et al. (2013), who assumed 
a absence of the dust torus. The unknown 7 -ray emission site 
may complicate the V 44 - b trend for FSRQs. Because of the 
larger photon energy density in the BLR, radiative cooling 
of electrons is stronger than in the dust torus. On the other 
hand, the relatively inefficient radiative cooling of electrons 
in the dust tours allows the electrons to be accelerated to 
higher energies by the inefficient stochastic mechanism. 

In our model, the radiation power is much greater than 
the magnetic-field and relativistic electron power, assuming 
that the baryon-loading (the ratio of the energy density of 
hadrons to electrons, 7 bi) is low, and will not have an im¬ 
pact on the total jet power {Pi+PB+Pe+Ph, where Ph is 
the power carried by hadrons) as long as rju < 10. Assum¬ 
ing r]hi = 1, we derive Ph = rjhiPe = 0.3 x lO'*® erg s“’^ and 
0.2 X10"^® erg s“^ for Period A and Period C, respectively. 3C 
279 has a black hole with mass (3-8) x 10 ®Mq (e.g., Gu et al. 
2001; Woo & 2002). The Eddington luminosity is therefore 
in the range (4-10) x 10“^® erg s“^. For low baryon loading, 
the total jet powers in our model are comfortably below the 
Eddington luminosity of 3C 279. 

Blazars are usually considered as cosmic-ray (GR) 
sources. Using the near-equipartition-log-parabola model, 
Dermer et al. (2015) gave the following expression for max¬ 
imum escaping proton energy 


Pmax(eV) = 1.4 x 10"°L®/^® ( A)1/8_ 

1^14 (( 

with /o = 1/3, fi = 10“^^"'®, and /2 


J-l/4 rl/8 

■^1 _ ( 6 ) 

1/4 *1/16 rll/16 ’ 1 1 

e Ss Jo 

= 10^/® . Using the 
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Figure 4. One-dimensional probability distribution of parameter values (left; the dashed curves are the mean likelihoods of samples and 
the solid curves are the marginalized probabilities) and two-dimensional contours of parameters (right; the regions enclosing 68%(95%) 
confidence level are shown) for Period C. The upper panel is for the input parameters, and the lower panel is for the output parameters. 


best-fit values for Period C (Table 1), we derive i?max ~ 
10^° — 2 X 10^° eV. However, it should be noted that 3C 
279 cannot contribute to the CRs with energy > 10^® eV 
at the Earth, because of the significant energy losses during 
the CRs traveling to us. 

Besides Emax, the other key quantity of a blazar is its 
CR luminosity Lcr. Secondary 7 rays are produced when 
high-energy CRs propagate towards us (e.g., Kalashev et al. 
2009, 2013), which depends primarily on the values of Lcr 
and Emax- Such kinds of non-variable secondary emission 
(variable on timescale of years; Prosekin et al. 2012) have 
been proposed to explain the steady VHE (> 100 GeV) 
emission from distant blazars (e.g., Essey et al. 2010, 2011; 
Murase et al. 2012; Aharonian et al. 2013; Takami et al. 
2013; Van et al. 2015). 

If 3C 279 emits high-energy cosmic rays, the secondary 
emission should be lower than its VHE emission, which 
puts an upper limit on Lcr. There are no recently re¬ 
ported VHE observations for 3C 279, so we use the sen¬ 
sitivity of MAGIC-H as the upper limit of secondary emis¬ 
sion. In Eig. 5, we show the secondary gamma-ray spectrum 
calculated by using the code TransportCR (Arisaka et al. 
2007; Gelmini et al. 2012; Kalashev & Kido 2014) and the 


EBL model of Franceschini et al. (2008) (see Finke et al. 
2010, for detailed comparisons for various EBL models). 
In the calculation, we assume a log-parabola GR distribu¬ 
tion with b — 1.48 and peak energy Epk = 5 x 10^® eV, 
and let the intergalactic magnetic field strength Bigmf = 
10“^® G (e.g., Essey, Ando & Kusenko 2011) and its coher¬ 
ence length Acoh = 1 Mpc. To make the secondary emis¬ 
sion lower than the sensitivity of MAGIG-II, the GR lumi¬ 
nosity is required to be < 6 x 10^® erg s“^, significantly 
lower than Lsyn and the apparent gamma-ray luminosity 
(~ lO'^® erg s“^). The upper limit of the corresponding ab¬ 
solute CR power Lj), is lO'^'^ erg s“^ (adopting = 28 
in Period C), a quarter of Pb {Pb is for a two-sided jet) in 
Period C. In Fig. 5, it can be seen that Cherenkov Telescope 
Array (CTA), having significantly improved sensitivity over 
the present imaging air Cherenkov arrays, will put stronger 
constraints on the CR luminosity. The upper limit assumes 
that ultra-high energy cosmic rays are not deflected and 
isotropized during transit, which is unlikely unless the blazar 
is found in structured regions, for example, clusters and fila¬ 
ments (Murase et al. 2012). Even if the UHECRs can freely 
escape from the source region, they can be deflected during 
transport across intergalactic space (Takami et al. 2014). 


MNRAS 000, 1-10 (2015) 































8 D. H. Yan, L. Zhang, S. N. Zhang 



'• Fermi data (Period Cjl 
MAGIC data 

r 

- MAGIC-II sensitivity 1 

r 

CTA sensitivity 

. secondary, 5E18 eV 7 

. 

■ • •. . 




v[Hz] 


Figure 5. Secondary 7 -ray emission produced in the propaga¬ 
tion of CRs. For comparison, the sensitivities of Imaging Atmo¬ 
spheric Cherenkov telescopes (MAGIC-II and CTA) and the non- 
simultaneous TeV MAGIC data for 3C 279 are shown. The sen¬ 
sitivities and TeV data are obtained through the SED Builder of 
ASDC (http://tools.asdc.asi.it). 
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2 orders of magnitude below that usually found in hadronic 
models (10-100 G) (Bottcher et al. 2013; Diltz et al. 2015). 

In this paper, we fitted the SEDs for Period A and Pe¬ 
riod C reported by Hayashida et al. (2015) with our one- 
zone leptonic model. This model can also explain the SED 
in the bright flare state of Period D with F{> 100 MeV) of 
10 “® photons cm“^ s“^ and tobs.var ~ 2 hours reported by 
Hayashida et al. (2015). In Fig. 6 , we show our best-fit result 
for the SED of Period D in Hayashida et al. (2015). In this 
fit, we use t 4 = 0.5 and Tdust = 800 K. The best-£t model re¬ 
quires ~ 10 and r ~ 0.2 pc. The deduced parameters are 
B ~ 0.8 G, (5d ~ 42, Ypk ~ 360 and i?' ~ 6 x 10^® cm (see the 
other parameters values in the caption in Fig. 6 ). The results 
of Period D, compared with those of Periods A and C, are 
consistent with b and 1/14 mutually increasing. In contrast, 
the HSP BL Lacs with large peak synchrotron frequencies 
have broader widths of the synchrotron SEDs (correspond¬ 
ing to smaller values of b) than FSRQs with smaller values 
of Us (Chen 2014). 

Because of the lack of X-ray spectrum and the poor data 
coverage of the synchrotron hump in Period B, we do not fit 
the SED in Period B. We found that the EG-dust emission 
in our model could account for the very hard 7 -ray spectrum 
(r.y ~ 1.7) in Period B with U 14 ~ 2.6, L 4 S 0.07, b 1.7, 
(/e ~ 3, t 4 0.5 and r ~ 0.3 pc; and we derived B ^ 0.3 G, 
(Id ~ 60, 7 pi; ~ 950, and Uduat 0.6 x 10“®erg cm“^. 

In conclusion, we have applied the MCMG technique 
developed for blazar studies by Yan et al. (2013) to the high- 
quality simultaneous data of 3G 279 (Hayashida et al. 2015), 
using a log-parabola EED (Dermer et al. 2014). We find that 
the 7 -ray SED is dominated by the dusty torus radiation, 
whereas the contribution of the BLR radiation field is very 
weak. This allows us to place the 7 -ray emission region out¬ 
side the BLR, at >0.1 pc, but within the IR radiation field 
of the torus. The reasonable fits obtained with curving par¬ 
ticle distributions are more simply expained with a stochas¬ 
tic acceleration mechanism. The allowed parameter range in 
3C 279 is well constrained at different epochs, and shows 
a positive correlation of the log-parabola width parameter b 
with peak synchrotron frequency Vs- Analysis of more FSRQ 
SEDs will be required to determine if this correlation is ro¬ 
bust. 


Figure 6. Best-fit model to the SED of 3C 279 (Hayashida et al. 
2015) for Period D. The parameter values are Uduat ~ 0.6 X 
10“^ erg cm~^, ~ 1.0, U 14 ~ 0.9, h ~ 1.4 and L 4 S ~ 0.3. 


If the hadrons in the jet are relativistic cosmic rays, as 
discussed above, these protons could, before escaping from 
the emission region, lose energies and produce secondary 
radiations via synchrotron and photohadronic interactions 
with low-energy photons (e.g., Dermer, Murase Inoue 
2014; Murase et al. 2014). However, the emission made by 
the protons in the emission region is negligible compared to 
the observed data, because the proton power (Ph) in our 
model is 2 - 3 orders of magnitude below the proton power 
required by hadronic model interpretations of FSRQ SEDs 
( 10 ^^- 10 ^® erg s“^) and the magnetic field in our model is 1 " 
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APPENDIX A: FITTING RESULTS WITH 
LARGER SIZE OF DUST TORUS 

We show the fitting results with a larger ruust, i-e., ruust = 
2 X 10^®(Ldisk/10^® erg cm (e.g., Ghisellini et al. 

2014). We then have ruust — 0.8 pc. Using ruust = 0.3, we 
derive Udust £ 2 x 10“"* erg cm“®. 

Figure Al shows the best-fit results with this larger ruust 
. One can see that the ht to the SED in Period A is com¬ 
parable to the fit showed in Section 3; however, the fit to 
lowest 7 -ray data in Period C is bad in the case of the larger 

■^dust- 

In Figs. A2 and A3, we show the one-dimensional prob¬ 
ability distribution of parameter values. The mean values 
and the marginalized 95% confidence intervals of the model 
parameters are reported in Table Al. In Period A, a larger 
(e and r is required in fitting with the larger rdust- The other 
input parameters change little, compared with those given in 
Section 3. The change in C,e leads to the significant changes 
in B' and (5 d. In Period C, the input parameters are nearly 
identical to those reported in Section 3. 

This paper has been typeset from a hU prepared by 

the author. 
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Figure Al. Best-fit models with a larger rdust to the SEDs of 3C 
279. Upper: for Period A; lower: for Period C. The inset shows 
the details of fit at X-ray energies. 
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Table Al. Input and output parameters values derived in the fittings with the larger r^ust- The mean values and the marginalized 95% 
confidence intervals (Cl) for interested parameters are reported. 


Input 


Ce 

b 

-^48 

*^14 

*4 

Cs 

r 

^dust 









(pc) 

(K) 

(lO^'^ erg S~-^) 

Period A 

3.25 

1.18 

0.19 

0.27 

5 

0.48 

0.7 

800 

0.15 

(95% Cl) 

1.52-4.91 

1.11-1.26 

0.17-0.20 

0.23-0.31 

- 

0.44-0.53 

0.1-1.2 


- 

Period C 

4.58 

1.50 

0.24 

0.52 

1 

0.79 

0.2 

800 

0.15 

(95% Cl) 

3.13-4.98 

1.43-1.56 

0.22-0.25 

0.48-0.57 

- 

0.71-0.87 

0.1-0.3 


- 

Output 


B 

■^D 

"^pk 

'^dust 

■^BLR 

R' 


Pt 



(G) 


(10~^ erg cm~^) 

(10~^ erg cm~^) 

(10^® cm) 

(lO**® erg s-1) 

(10“^® erg s“^) 


Period A 

0.5 

23 

290 

0.1 

0.005 

3 

0.9 

0.6 


(95% Cl) 

0.3-0.8 

19-26 

260-320 

0.04- 0.25 

0.001-0.3 

- 

- 

- 


Period C 

1.0 

29 

300 

0.21 

0.6 

0.8 

0.5 

1.0 


(95% Cl) 

0.9-1.3 

27-31 

290-310 

0.20- 0.22 

0.07-2 

- 

- 

- 




Figure A2. One-dimensional probability distribution of param¬ 
eter values for Period A. 


Figure A3. One-dimensional probability distribution of param¬ 
eter values for Period C. 
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